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We study the quasi-elastic contribution to the inclusive {vi,viN), {vi,l~N), (ui,i'iN) and 
{vi,l'^N) reactions in nuclei using a Monte Carlo simulation method to account for the rescat- 
tering of the outgoing nucleon. As input, we take the reaction probability from the microscopical 
many body framework developed in Phys. Rev. CTO (2004) 055503 for charged-current induced 
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^ ■ I. INTRODUCTION 

in : 

' Neutrino physics is at the forefront of current theoretical and experimental research in astro, nuclear, and particle 
^ . physics. Indeed, neutrino interactions offer unique opportunities for exploring fundamental questions in these domains 
' of the physics. One of these questions is the neutrino-oscillation phenomenon, for which there have been conclusive 
, positive signals in the last years [1]. Neutrino-oscillation experiments are presently evolving from the discovery to 
fSj the precision phase. This new generation of experiments faces a major difficulty: the elusive nature of the neutrinos. 
[ The presence of neutrinos, being chargeless particles, can only be inferred by detecting the secondary particles they 
• create when colliding and interacting with matter. Nuclei are often used as neutrino detectors, thus a trustable 
1^ interpretation of neutrino data heavily relies on detailed and quantitative knowledge of the features of the neutrino- 
I nucleus interaction [2]. For instance, in the case of neutrino processes driven by the electroweak Neutral Current (NC), 
' the energy spectrum and angular distribution of the ejected nucleons are the unique observables. There is a general 
consensus among the theorists that a simple Fermi gas model, widely used in the analysis of neutrino oscillation 
O • experiments, fails to provide a satisfactory description of the measured cross sections, and inclusion of further nuclear 
<]J I effects is needed [3]. 

' Simultaneously, in recent years there have also been some initiatives aiming at understanding the quark and gluon 
' substructure of the nucleon. The flavor dependence of the nucleon's spin is a significant fundamental question that 
is not yet fully understood. Experiments [4, 5] measuring the spin content of the nucleon have reported conflicting 
results on the amount of nucleon spin carried by strange quarks [6] . Recently, the FINeSSE collaboration at Fermilab 
has suggested [7, 8] that Quasi-Elastic (QE) neutrino-nucleus scattering, observed using a novel detection technique, 
provides a theoretically clean measure of this quantity. In this context, it is also necessary to control nuclear effects. 

At intermediate energies, above the nuclear giant resonance and below the A(1232) regions^, neutrino-nucleus 
interactions have been studied within several approaches. Several different Fermi gas. Random Phase Approximation 
(RPA), shell model and superscaling based calculations have been developed during the last 15 years [19]-[38]. Some 
of these approaches have been also employed to compute neutrino or antineutrino induced single-nucleon emission 
cross sections. Most of the calculations use the plane wave and distorted wave impulse approximations (PWIA and 
DWIA, respectively), including or not relativistic effects. The PWIA calculations neglect all types of interactions 
between the ejected nucleon and the residual nuclear system, and therefore such a framework constitutes a poor 
approximation to evaluate nucleon emission cross sections. However, the PWIA has been often used to compute the 
ratio of proton (f, p) to neutron (i^, n) yields, which at low neutrino energies and for light nuclei might be rather 
insensitive to rescattering effects. 



X 



There exists an abundant literature studying these two regions. See for instance a recent paper [9] on the excitation of nuclear giant 
resonances in neutrino scattering off nuclei or the older works of Refs. [10]-[15] also studying QE neutrino-nucleus scattering at low 
energies. On the A— excitation in neutrino reactions, there are works, among others, by Alvarez-Ruso and collaborators [16], S.K. Singh 
and collaborators [17] and Lalakulich y Paschos [18] . 
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Within the DWIA the ejected nuclcon is described as a solution to the Dirac or Schrodinger equation with an optical 
potential obtained by fitting elastic proton-nucleus scattering data. The imaginary part accounts for the absorption 
into unobserved channels^. This scheme, first developed in (e, e'p) studies where the final nucleus is left in the ground 
or in a particular excited state, is incorrect to study nucleon emission processes where the state of the final nucleus 
is totally unobserved, and thus all final nuclear configurations, either in the discrete or on the continuum, contribute. 
The distortion of the nucleon wave function by a complex optical potential removes all events where the nucleons 
collide with other nucleons. Thus, in DWIA calculations, the nucleons that interact are lost when in the physical 
process they simply come off the nucleus with a different energy, angle, and maybe charge, and they should definitely 
be taken into account. A clear example which illustrates the deficiencies of the DWIA models is the neutron emission 
process: {vi,l~n). Within the impulse approximation neutrinos only interact via Charged Current (CC) interactions 
with neutrons and would emit protons, and therefore the DWIA will predict zero cross sections for CC one neutron 
knock-out reactions. However, the primary protons interact strongly with the medium and collide with other nucleons 
which are also ejected. As a consequence there is a reduction of the flux of high energy protons but a large number 
of secondary nucleons, many of them neutrons, of lower energies appear. 

The distortion by a real potential does not eliminate the events where there are nucleon collisions. But, it does not 
account either for the changes of energy, direction and charge of the nucleons induced by these collisions. Besides, it 
cannot account either for more than one nucleon knock-out events. 

In this work, we study the QE contribution to the inclusive {vi,viN), {vi,l~N), {i>i,PiN) and (PiJ^N) reactions 
in nuclei. We use a Monte Carlo (MC) simulation method to account for the rescattering of the outgoing nucleon. A 
reliable description of the gauge bosons (W^ and ) absorpt ion in the nucleus is the first essential ingredient. 

For CC driven processes, we use the many body framework developed in Ref. [34]. Starting from a Local Fermi 
Gas (LFG) picture of the nucleus, which automatically accounts for Pauli blocking, several nuclear effects are taken 
into account in the scheme of Ref. [34]: i) a. correct energy balance, using the experimental values, is enforced, 
ii) Coulomb distortion of the charged leptons is implemented by using the so called "modified effective momentum 
approximation" [40], in) medium polarization (RPA), including A— hole degrees of freedom and explicit pion and rho 
exchanges in the vector-isovector channel of the effective nucleon-nucleon force, and Short Range Correlation (SRC) 
effects are computed, and finally iv) the nucleon propagators are dressed in the nuclear medium, which amounts to 
work with a LFG of interacting nucleons and it also accounts for reaction mechanisms where the gauge boson, 
or W~, is absorbed by two nucleons (the real part of the nucleon selfenergy modifies the free nucleon dispersion 
relation, while the imaginary part takes into account two nucleon absorption reaction channels). The model has no 
free parameters. The W^N couplings and form factors are fixed in the vacimm, while the main features concerning the 
nuclear corrections, expansion parameter and all sorts of constants, are completely fixed from previous hadron-nucleus 
studies (pionic atoms, elastic and inelastic pion-nucleus reactions, A— hypernuclei, etc.) [41-43]. Thus, the model is a 
natural extension of previous studies [41 45] on electron, photon and pion dynamics in nuclei, and should be able to 
describe inclusive CC QE neutrino and antineutrino nuclear reactions at intermediate energies of interest for future 
neutrino oscillation experiments. Even though the scarce existing CC data involve very low nuclear excitation energies, 
for which specific details of the nuclear structure might play an important role, the model of Ref. [34] provides one 
of the best existing combined description of the inclusive muon capture in ^^C and of the measurements of the ^^C 
{h>^,,^ir)X and ^^C {ve,e^)X reactions near threshold. Inclusive muon capture from other nuclei is also successfully 
described by the model. Besides, above, let us say 80 or 100 MeV of energy transferred to the nucleus, this many 
body framework leads also to excellent results for the (e, e') inclusive reaction in nuclei, not only in the QE region, 
but also when it is extended to the study of the A— peak and the dip region (situated between the QE and the A 
peaks) [44]'^ and to the description of the absorption of real photons by nuclei [45]. 

In Sect. II, we extend the model of Ref. [34] to NC driven processes, both for neutrino and antineutrino induced 
nuclear reactions in the QE region. Thus we compute, for a fixed incoming neutrino or antineutrino Laboratory (LAB) 
energy, the inclusive QE cross sections (f a / dft' dE' ($!', E' arc the solid angle and energy of the outgoing lepton) for 
(yi^vi), (i>i,l~), {vi,Di) and (i/;,Z+) processes. This cross section gives us the reaction probability and it is the first 
required ingredient to start with our cascade model to describe the collisions suffered by the nucleons through their way 
out of the nucleus'*. Details on the MC simulation are given in Sect. III. The MC method used here was designed for 



^ For nuclcon energies above 1 GcV, the Glauber model [39] , which is a multiple-scattering extension of the eikonal approximation, has 
also been used (sec for instance the recent work of Ref. [35]). In this approach, a relativistic plane wave is modulated by a factor which 
accounts for the absorption into unobserved channels. 

^ Data in ^^C, '"'Ca and ^OSpj^ differential cross sections for different electron kinematics and split into longitudinal and transverse 
response functions are successfully described. 

Besides, we also compute differential cross sections with respect to cPr. Thus, we also know the point of the nucleus where the gauge 
boson was absorbed, and we can start from there our MC propagation of the ejected nucleon. 
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single and multiple nucleon and pion emission reactions induced by pions [42, 46] and has been successfully employed 
to describe inclusive (7,7r), (7,7V), (7, iViV),..., {■y,NTT),... [47, 48], (e,e'7r), (e,e'iV), (e, e'A^A^),..., (e, c'A^tt),... [49] 
reactions in nuclei or the neutron and proton spectra from the decay of A hypernuclei [50]. Thus, we are using a quite 
robust and well tested MC simulator. 

In Sect. IV we discuss our results and the main conclusions of this work are outlined in Sect. V. We start presenting 
(Subsect. IV A) results for the inclusive QE NC cross sections for both neutrino and antineutrino beams in several 
nuclei and in the next subsection we show results for inclusive {i/i,i/iN), [vi, l~N), {i>i,i/iN) and {i>i,l~^N) reactions in 
nuclei at low energies, obtained from our cascade model. Some preliminary results of this work were already presented 
in Rcf. [51]. Finally, in the appendix we give explicit expressions for the NC nucleon tensor, both in the impulse 
approximation and when RPA corrections are taken into account. 

To end this introduction, we would like to devote a few words on the applicability of the nuclear model used here. 
One might think that a LFG description of the nucleus is poor, and that a proper finite nuclei treatment is necessary. 
For inclusive processes and nuclear excitation energies of around 100 MeV or higher, the findings of Refs. [43], [44] 
and [45] clearly contradict this conclusion. The reason is that in these circumstances one should sum up over several 
nuclear configurations, both in the discrete and in the continuum, and this inclusive sum is almost not sensitive to 
the details of the nuclear wave function^ , in sharp contrast to what happens in the case of exclusive processes where 
the final nucleus is left in a determined nuclear level. On the other hand, the LFG description of the nucleus allows 
for an accurate treatment of the dynamics of the elementary processes (interaction of gauge bosons with nucleons, 
nucleon resonances, and mesons, interaction between nucleons or between mesons and nucleons, etc.) which occur 
inside the nuclear medium. Within a finite nuclei scenario, such a treatment becomes hard to implement, and often 
the dynamics is simplified in order to deal with more elaborated nuclear wave functions. This simplification of the 
dynamics cannot lead to a good description of nuclear inclusive electroweak processes at the intermediate energies of 
interest for future neutrino experiments. 



II. EXTENSION OF THE MODEL OF REF. [34] TO NEUTRAL-CURRENTS 

A. General formalism: hadronic tensor and many body expansion 

We will first focus on the neutrino induced inclusive reaction driven by the electroweak NC 

ui{k) + Az ui{k') + X (1) 

and we will follow the same notation and convention as in Ref . [34] . The double differential cross section, with respect 
to the outgoing neutrino kinematical variables, for the process of Eq. (1) and for massless neutrinos, is given in the 
LAB frame by 

' ' ' 2Wi sin^ — + W2 cos^ — - W3 ' ' ' ' ' sin^ - } (2) 




dfl{k')d\k'\ 47r2 1 2 2 M, 



with k and k' the LAB neutrino momenta, G = 1.1664 x 10~^^ MeV"^, the Fermi constant, 9' the outgoing neutrino 
scattering angle and M, the target nucleus mass. To obtain Eq. (2) we have neglected the four-momentum carried 
out by the intermediate .Z— boson with respect to its mass. Finally, the three independent, Lorentz scalar and real, 
structure functions, Wi{q'^), enter into the definition of the hadronic tensor, W'^'^, which includes all sort of non- 
leptonic vertices and corresponds to the neutral current electroweak transitions of the target nucleus, i, to all possible 



^ The results of Ref. [34] for the inclusive muon capture in nuclei through the whole periodic table, where the capture widths vary from 
about 4x10'' s-i in to 1300 xlO* s'^ in '^°^Ph, and of the LSND measurements of the ^^C (^'p,/^-)X and ^^C (i^e,e-)X reactions 
near threshold indicate that the predictions of our scheme, for totally integrated inclusive observables, could even be extended to much 
smaller, of the order of 10 or 20 MeV, nuclear excitation energies. In this respect, the works of Refs. [52] and [53] for inclusive muon 
capture and radiative pion capture in nuclei, respectively, turn out to be quite enlightening. In those works, continuum shell model 
results are compared to those obtained from a LFG model for several nuclei from ^^C to ^°*Pb. The differential decay width shapes 
predicted for the two set of models are substantially different. Shell model distributions present discrete contributions and in the 
continuum appear sharp scattering resonances. Despite the fact that those distinctive features do not appear in the LFG differential 
decay widths, the totally integrated widths (inclusive observable) obtained from both descriptions of the process do not differ in more 
than 5 or 10%. The typical nuclear excitation energies in muon and radiative pion capture in nuclei are small, of the order of 20 MeV, 
and thus one expects that at higher excitation energies, where one should sum up over a larger number of nuclear final states, the LFG 
predictions for inclusive observables would become even more reliable. 



4 



final states. It is thus given by (in our convention, we take eoi23 = +1 and the metric g^^"^ = (+, — , — , — ))^ 

^rflo■ I 



2Mi AMf 



^(27r)-^54(p; -p- qmiiAomifbnoimr 



f 



^ ^ M? 2M2 '^'^3+ M2 2M2 2M2 ^ 

with P^^ the four-momentum of the initial nucleus {Alf = P^), Pj- the total four momentum of the hadronic state / 
and q = k — k' the four momentum transferred to the nucleus. The bar over the sum denotes the average over initial 
spins. By construction, the hadronic tensor accomplishes 

Wf" = W^" + iW^" (4) 

with W^"^ i^a'^) I'eal symmetric (antisymmetric) tensors, and finally for the NC we take 

i^, = ^^T^Cl - I sin^ ew - 75)*« - *d7''(l + ^ sin^ - 75)*d - "^s^^il + ^ sin^ Ow - 75)*« (5) 

with ^frf and \E's quark fields, and 9w the Weinberg angle (sin^ 9w = 0.231). 

Taking q'm the z direction, ie, q = \q\uz, and P^^ = (A/^, 0), it is straightforward to find the six structure functions in 
terms of the W°°, W^^ = W^^ , W^^, W^^ and W^'' components of the hadronic tensor^. The neutrino cross section, 
Eq. (2), does not depend on Mj, as can be seen from the relations of Eq. (6), and also note that the structure functions 
^4,5,6 do not contribute in the limit of massless neutrinos. 

The cross section for the antineutrino induced nuclear reaction 

n{k) + Az ^ vi{k') + X (7) 

is easily obtained from Eq. (2), just by changing the sign of the parity- violating term {W-^). 

The hadronic tensor is determined by the Z"— boson sclfenergy, n^''(g), in the nuclear mcdimn. We follow here the 
formalism of Ref. [34], and we evaluate the selfenergy, S^(A:; p), of a neutrino, with four-momentum k and helicity r, 
moving in infinite nuclear matter of constant density p. We find, 

wr = -0(9°) (^^^^ 7 i"^'' + "^'^ f^'^^ 

wr = -e(,°) 7 ^ f""" ~ """^ ^"'^ '^'^^^ 

where we have used the Local Density Approximation (LDA) to obtain results in a finite nucleus of density p{r), and 
g is the gauge weak coupling constant, g = e/ sin 6w, related to the Fermi constant: Gjsf^ = g'^/SMy^, with e the 
electron charge. 

As we see, the basic object is the selfenergy of the Gauge Boson (Z^) inside of the nuclear medium. As it is done in 
Ref. [44] for electro induced nuclear reactions, we plan to perform a many body expansion, where the relevant gauge 
boson absorption modes would be systematically incorporated: absorption by one nucleon, or a pair of nucleons or 
even three nucleon mechanisms, real and virtual meson (tt, p, • • •) production, excitation of A of higher resonance 
degrees of freedom, etc. In addition, nuclear effects such as RPA or SRC should also be taken into account. Some of 
the absorption modes are depicted in Fig. 1. 



Note that: (i) Eq. (3) holds with states normaUzcd so that {p]p') = {2iT)''"2p[)5'^ (p — p'), (ii) the sum over final states / includes an 
integration J (2^)^^2g ' ^"'^ each particle j making up the system /, as well as a sum over all spins involved. 
Thus, one readily finds 

^1 = 1771-. ^^2 = — VK»"-^VV'"" + -^(W'"-H/-")-2-^ReTyOM , Ws = -i—- (6) 
2Mi 2Mi \ \q\^ \q\ J I? I 
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FIG. 1: Diagrammatic representation of some diagrams contributing to the Z— self energy. 



B. QE contribution and Pauli blocking 

The virtual Z° can be absorbed by one nucleon leading to the QE contribution of the nuclear response function. 
Such a contribution corresponds to a Iplh nuclear excitation (first of the diagrams depicted in Fig. 1). To evaluate 
this selfenergy, the free nucleon propagator in the medium is required. 

S{p ■,p) = i^ + M)G{p ; p), G{p ; p) = _ ^ + ^^^iP° " E{pmkF - \p |)) (10) 

with the local Fermi momentum kpir) = {Sn"^ p{r)/2y/^ , M = 940 MeV the nucleon mass, and E{p) = sjM"^ 
Wc will work on a non-symmetric nuclear matter with different Fermi sea levels for protons, A;^, than for neutrons, 
A:^ (equation above, but replacing p/2 by p^ or p„, with p — p^^ On the other hand, for the Z^NN vertex we 
take (-/V = n or p) 

<N;p'=p + q\j^M\N;p >= u(p '){VS - A%)uip) (11) 
with spinor normalization given by uu = 2m, and vector and axial nucleon currents given by 

y;^ = 2x (Ff(g2)7«+vz^^§pc7-g.)^, A%={Gi{q')r7^ + GUq')q''l5)^ (12) 

Invariancc under G-parity has been assumed to discard a term of the form (p^ + p'^)75 in the axial sector, while 
invariance under time reversal guarantees that all form factors are real. Besides, and thanks to SU(3) symmetry, 
some relations exist among the NC form factors and the CC Ga and Gp, see Ref. [34]) and the electromagnetic 
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ones^ (J'f'", HpFl^ and ^inF^) 

{FfY'"' = ±FY - 2 sin^ e^Ff - (13) 

{^izF^Y''" = ±^ivF^ -2^in^ewl^p,nFr -\^J^sF^ (14) 

(GIp)'''" = ±G^,p-G1p (15) 

where Fi,iisF2,G\ and Gp are the strange vector and axial nucleon form factors [55]. The pseudoscalar part of the 
axial current does not contribute to the differential cross section for massless neutrinos and for the rest of strange 
form factors we use the results of the fit II of Ref . [56] , 

with Qj^ = —0.15 and M\ = 1049 MeV. With all of these ingredients is straightforward to evaluate the contribution 
to the selfenergy of the first diagram of Fig. 1, which leads to 

X E {Q{k^{r)-\pm\p + q\-k^{r))A''^^{p,q)] Siq° - Q + E{p) - E{p + q)) (17) 

where Q is the experimental Q— value, included in order to properly reproduce the energy threshold. For inclusive 
observables we have set it to zero to approximately take into account the possibility of elastic scattering and transitions 
to nuclear excited states. The (fip integrations above can be analytically done and all of them are determined by the 
imaginary part of the relativistic Lindhard function, U R{q,kp ,kp). Explicit expressions can be found in Appendix 
B of Ref. [34]. The NC nucleon tensor, A^, can be found in Appendix Al. The non-relativistic reduction of the 
hadronic tensor can be obtained by replacing the factors M/E{p) and M/E{p + q) in Eq. (17) by one. Explicit 
expressions can be found in Appendix C of Ref. [34] . 

To finish this section, we devote a few words to the Low Density Theorem (LDT). At low nuclear densities the 
imaginary part of the relativistic Lindhard function can be approximated by 

ImUniq, , ) ~ -^PN-^S{q° + M- E{q)) (18) 

and thus one readily finds 

<^iyi+Az^iyi+X ~ Nffi^i+n^^l+n + Zai,i + p^i,l+p, N = A — Z (19) 

which agrees with the LDT. 

C. QE contribution: RPA nuclear correlation and FSI effects 

Pauli blocking, through the imaginary part of the Lindhard function, is the main nuclear effect included in the 
hadronic tensor of Eq. (17). RPA nuclear correlation and FSI effects played a crucial role for inclusive QE CC 
neutrino-nucleus reactions, and we include those here, following the same formalism as in Ref. [34]. 

1. RPA 

We replace the Iplh contribution to the selfenergy by an RPA response as shown diagrammatically in Fig. 2. 
For that purpose wc use an effective ph-ph interaction of the Landau-Migdal type, with explicit contribution of pion 
and rho meson exchanges, 

V = Co|/o(p) + /o(p)nr2 + go{p)aia2^ + TiT2^Vi{q)aiqa2q + Vt{q){aia2 - aiqa2q)^ (20) 



We use the parameterization of Galster and collaborators [54], which is also compiled in Ref. [34]. 
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where (T and r are Pauli matrices acting on the nucleon spin and isospin spaces, respectively, and q = q/\q\- We take 
the coefficients /o, /q and from Ref. [57], 



fMr)) = ^/r^ 



+ 



1- 



p{t) 
P(0) 



(ex) 



(21) 



where 



(in) 



Jo 

f(ex) _ 

Jo — 

(in) 

96 



0.07 



-2.15 

5^'"' = 90 = 0.575 



fl(ex) 

Jo 

fl(in) 

Jo 



: 0.45 

0.33 



(22) 



and Co = 380 MeVfm"^. In the S = 1 = T channel (ctcttt operator) we use an interaction with explicit tt (longitudinal) 
and p (transverse) exchanges, which has been used for the renormalization of the pionic and pion related channels in 
different nuclear reactions at intermediate energies [42]-[45] . The strengths of the ph-ph interaction in the longitudinal 

and transverse channel are given by 



mi 



9'M) 



P 

— = 0.08, 

47r 



= 1200 MeV 




V 



/ V 
Z ( 

FIG. 2: Set of irreducible diagrams responsible for the polarization (RPA) effects in the Iplh contribution to the selfenergy. 
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yt{q\ «1 = ^ jC-p J + j ' = 2, = 2500 MeV, m, = 770 MeV (23) 

The SRC functions g[ and g[ have a smooth q— dependence [41, 58], which we wiU not consider here^, and thus we 
will take g'i{q) = g'tiq) = g' = 0.63 as it was done in the study of inclusive nuclear electron scattering carried out in 
Ref. [44], and also in some of the works of Ref. [42, 43]. 

The above interaction has been successfully tested in different nuclear processes at intermediate energies [41]-[49], 
and we will use here the same form and parameters as in our recent study of CC neutrino nucleus reactions [34]. 
As it is explained there, A(1232) degrees of freedom are also taken into account. Given the spin-isospin quantum 
numbers of the A resonance, these degrees of freedom only modify the vector-isovector [S = 1 = T) channel of the 
RPA response function. The ph-Ah and Ah-Ah effective interactions are obtained from the interaction of Eq. (20) 
by replacing — > S, r — > T, where S, T are the spin, isospin NIS. transition operators [41] and / f* = 2.13 /, for 
any A which replaces a micleon. 

Thus, the V lines in Fig. 2 stand for the effective ph(Ah)-ph(Ah) interaction described so far. Keeping track of the 
operators responsible, we have examined and renormalized all different contributions to the NC nucleon tensor A'^ , 
by summing up the RPA series depicted in Fig. 2. The procedure is discussed in detail in Ref. [34], now however the 
isospin structure of the Z^NN vertex 



with r'' some matrix in the Dirac space and F some form factor, does not select only the isovector channels of the in 
medium effective ph(Ah)-ph(Ah) interaction, as it is the case for the CC induced processes studied in [34], and hence 
we also find contributions from the isoscalar part of the effective force. 

For simplicity and to compute these polarization corrections, we have assumed a symmetric nuclear matter with 
the same Fermi sea level for protons and neutrons. On the other hand, since the ph(Ah)-ph(Ah) effective interaction 
is non-relativistic, we have computed polarization effects only for the leading and next-to-leading terms in the p/M 
expansion. Thus, order O (^kFP^/M'^,kFp''^/M^,kpq°/M'j has been neglected. We have made an exception to the 
above rule, and since /Uz could be relatively large, we have taken fizF2\q\/M to be of order 0(0) in the p/M 
expansion. With all these ingredients, we find SW^^pj^, contribution which has to be added to the hadronic tensor, 
W^"^, given in Eq. (17) to account for the medium polarization effects, 

= -5i^r*'-1^e(,»)/||,3|L_J£_,_.),(,o^^ 

X {e{kf{r)-\pne{\p + q\-kArnSA'S,^{p,g)} (25) 

L J p"=E{j>) 

The 00, Oz, zz, XX and xy components of the RPA contribution to the NC nucleon tensor, 5A'^j^{p,q), are given in 
Sect. A 2 of the Appendix. There, the density dependent tensor SA'^j^{p, q) is given in terms of the parameters of the 
effective ph(Ah)-ph(Ah) force and U{q,kp) = Un + Ua, the Lindhard function of Ref. [58], which for simplicity we 
evaluate^" in an isospin symmetric nuclear medium of density p. The different couplings for N and A are incorporated 
in Un and Ua and then the same interaction strengths Vi and Vt are used for ph and Ah excitations [41]-[49]. Note 
that, backward (crossed term of the Lindhard function) propagating ph and Ah excitations are also taken into account 
within our framework. For positive values of the backward propagating ph excitation has no imaginary part, and 
for QE kinematics Ua is also real. 

Finally, we should stress that the /g, /o and go terms of the effective interaction cannot produce Ah excitations 
and therefore, when these terms are involved in the RPA renormalization, only the nucleon Lindhard function (C/jv) 
appears (see Eqs. (A9)). 



2. FSI: dressed nucleon propagators in the nuclear medium 

Once a ph excitation is produced by the virtual Z*^— boson, the nucleon is interacting with the rest of nucleons of 
the nucleus, colliding many times, thus inducing the emission of other nucleons. The result of it is a quenching of the 



This is justified because taking into account the dependence leads to minor changes for low and intermediate energies and momenta, 
where this effective ph-ph interaction should be used. 

The functions Uff and C/a are defined in Eqs. (2.9) and (3.4) of Ref. [58], respectively. 
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QE peak respect to the simple ph excitation calculation and a spreading of the strength, or widening of the peak. 
In our many body scheme we will account for the FSI by using nucleon propagators properly dressed with a realistic 
selfenergy in the medium, which depends explicitly on the energy and the momentum [59]. Hence we substitute 
the particle nucleon propagator, G{p;p), in Eq. (10) by a renormalized nucleon propagator, Gfsi(p;p)) including the 
nucleon selfenergy in the medium, E(p°,p ; p), 

Gfsi(p;p) = -5 ^, ^ . (26) 

with E{p) = M + p"^ /2M. This approach led to excellent results in the study of inclusive electron and CC neutrino 
scattering from nuclei [34, 44]. Since the model of Ref. [59] is not Lorentz relativistic and it also considers an 
isospin symmetric nuclear medium, we will only discuss the FSI effects for nuclei with approximately equal number of 
protons and neutrons, and using non-relativistic kinematics for the nucleons. Thus, we have obtained Eq. (26) from 
the non-relativistic reduction of G(p; p), in Eq. (10), by including the nucleon selfenergy. 

To account for FSI effects in an isospin symmetric nuclear medium of density p we should make the following 
substitution [34, 44] 



2e(g") j j^eikAr) - \P m\p + q 1 - kF{rm-7T)S{q° + E{p) - E{p + q))A-'^{p, <z)l,o=^(,-) 

e(g°) 

47r2 



/ d'p r ckjSh{u;,p; p)5p(90 + u;,p + q; pJA'^^'ip, q\o^Eif) , = ^ ^<pa (27) 

in the expression of the hadronic tensor (Eqs. (17) and (25) ). Sh,Sp are the hole and particle spectral functions 
related to nucleon selfenergy E by means of 

Sp,h{'^,p;p) = : r2 — T2 (28) 

a; — S(p) — ReS(a;,p; p) + ImE(a;,p;p) 

with a; > or w < /X for 5p and Sh, respectively. The chemical potential p, is determined by 

M = M+^ + ReS(M,fc^) (29) 

The d^p integrations have to be done now numerically, and since the imaginary part of the nucleon selfenergy for 
the hole states is much smaller than that of the particle states at intermediate nuclear excitation energies, we make 
the approximation of setting to zero ImS for the hole states. Thus, we take 

Sh{ioJ; p) = 5{u - E{p))eip - E{p)) (30) 

where E(j)) is the energy associated to a momentum p obtained self consistently by means of the equation 

E{p) = E{p) + ReS(^(2?),p;p) (31) 

The same approximation is also used in some of the calculations for the particle spectral function. The effects of this 
approximation will be discussed in Sect. IV B. 



III. THE MONTE CARLO SIMULATION 



A. Kinematics of the outgoing nucleon in the first step 



In a previous work [34] and in Sect. II we carried out a thorough evaluation of the CC and NC inclusive neutrino 
and antineutrino induced nuclear reactions in the QE region. Thus we have determined, for a fixed incoming neutrino 
or antineutrino LAB energy, the inclusive QE cross section (P a / dfl' dE' (il', E' are the solid angle and energy of the 
outgoing lepton). The absorption of the gauge boson {W^ or Z^) with four momentum by one nucleon constitutes 
the reaction mechanism (Iplh excitation) and the corresponding reaction probability is determined by d^a/dfl'dE' . 



This is not entirely eorrect, sinee two nueleon absorption modes are also considered, when the FSI effects described in Sect.IIC2 are 
taken into account. We will discuss this point later in Subsect. IV B. 
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Moreover, within our scheme, we obtain (Pa/dft'dE' after performing an integration over the whole miclear volume 
(see Eqs. (16) and (17) of Ref. [34] and Sect. II, respectively). Thus, for a fixed transferred four momentum q^, chosen 
according to cfa/dil'dE' , we can randomly select the point of the nucleus where the absorption takes place using the 
profile d^a/dn'dE'd^r evaluated in Ref. [34] for CC and in Sect. II for NC. 

Now, we need on top of that to have the distribution of three-momenta of the outgoing nucleon. This can be done 
by not performing the integration over the three momentum, p, of the nucleon occupied states in [34] and in Sect. II 
(see for instance Eq. (17) for the the NC case), since the outgoing nucleon momentum is p + q. However, in order 
to incorporate the nucleon scattering with other nucleons, we choose the alternative procedure of generating events 
probabilistically one by one, with a weight given by d^a/dfl'dE'd^r. For each event in a certain position rand with a 
transferred four momentum we generate a random momentum p from the local Fermi sea. The vector p + q gives us 
the direction of the nucleon. The energy of the nucleon is then obtained by imposing energy conservation assuming in 
this step and throughout all the MC simulation that the nucleons move in a Fermi gas under an attractive potential 
equal to the local Fermi energy and therefore 

E' = E{p) + q'' (32) 

where E{p) = \/]P + M^ — kp{r)/2M . This will provides us with the modulus of the outgoing momcntum-'^^ p' . If it 
happens that \p'\ < kpir) (the local Fermi momentum) then the event is Pauli blocked, it is dismissed and another 
event is generated. Thus, we have already the configuration of the final state after the first step, namely, one nucleon 
produced in the point r of the nucleus with momentum p' . With respect to having a proton or a neutron in the final 
state, this is trivially done: for the CC case, the outgoing nucleon is a proton (neutron) for neutrino (antineutrino) 
induced processes, while for NC, in Sect. II, the reaction probability was already split into a proton and a neutron 
induced ones (see for instance Eq. (17))^^. 



B. Nucleon propagation 

The nucleons in the nucleus move under the influence of a complex optical potential. The imaginary part of the 
potential is related to the probability of nucleon quasielastic collisions in the nucleus (and extra pion production 
at higher energies, which we do not consider here). We consider explicitly these collisions since they generate new 
nucleons going outside the nucleus. As with respect to the real part, we use it to determine the classical trajectories 
that the nucleons follow in the nucleus between collisions. 

As done in Refs. [48, 49] we take as the real part of the nucleon-nucleus potential 

Vir) = - Sir) = -§ = -^ {r'Pi^^f' (33) 

It represents the interaction of a single nucleon with the average potential due to the rest of the nucleons. This choice 
of V{r) means that the total nucleon energy is the difference between its kinetic energy and the Fermi energy, kp/2M . 

With respect to collisions in our MC simulation we follow each excited nucleon by letting it move a short distance 
d such that Pd « 1 {P represents the probability per unit of length for a quasielastic collision ). The new position 
(r ') and momentum (p') are taken from the Hamiltonian equations as 

r' = f+Sf=f+-^d 

\P\ 

^1 ^ dVE(p)d f 

P = P + Sp; Sp=-- ^—jz^ 34 

or p \r\ 

which follow from the Hamilton equations or equivalently from energy and angular momentum conservation. 

Our code selects randomly, according to the reaction probabilities which will be discussed in Subsect. IIIC, if the 
nucleon is scattered or not and, in the case of scattering, what kind of process takes place. If no collision takes 
place, we move the nucleon again. When the nucleon leaves the nucleus we stop the process and it is counted as a 
contribution to the cross section. If a NN scattering is selected instead, we take a random nucleon from the Fermi 



Indeed, we have = (E' + A:|,(r)/2M)2 - M^. 

Note that polarization corrections (Eq. (25)) were estimated in isospin symmetric nuclear matter, and thus approximately the RPA 
term equally contributes to both the proton and neutron Z'^ absorption channels. 
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sea and calculate the initial kinematical variables (P^ and s, full four- momentum of the nuclcon-nucleon system in 
the nuclear frame and invariant energy, respectively). Then, a cos is selected, according to the expression given 
in Eq. (37) below and taken from Ref. [48]. This expression gives us the correct probability given by da^^ / (Klc.m. 
plus Pauli blocking restrictions. We take also into account Fermi motion and renormalization effects in the angular 
dependence. We take them into account by multiplying each event by a weight factor 

? = A(7Vi,7V2)a^i^='pjv. (35) 

where l/\{Ni,N2) is the probability by unit of length that one nucleon iVi collides with another nucleon N2. In 
average, this factor ^ is equal to one. Explicit formulas are given in the next subsection. 

Our method assumes that the nucleons propagate semiclassically in the nucleus. The justification of this hypothesis 
for reactions induced by real photons is given in Ref. [48], and has also been successfully used in Refs. [46, 50]. In the 
next section we give some detail on the evaluation of the equivalent NN cross section in the medium. 



C. NN cross sections 



We are using the parameterization of the NN elastic cross section given in the Appendix of Ref. [48]. Since for 
particles of low momenta, the MC induces large errors, we are not considering collisions of nucleons with kinetic 
energies below 30 McV. This is to say, we do not follow the path through the nuclear medium of nucleons with kinetic 
energies below 30 MeV, and we just consider that those nucleons get out of the nucleus without suffering further 
collisions, in which eventually they could change charge or loss some more energy. 

On the other hand, the reaction probability will change due to the nuclear medium effects (Fermi motion, Pauli 
blocking and medium renormalization). Then, according to Ref. [48], the expression for the mean free path (A) of the 
nucleon is given by 



1 



\{Ni) 



Q{kl{r)-\p2\) 



{A - Z) 



.Nm 



is) 



\p 



llab 



1^1 



(36) 



with P^, the full four-momentum of the NN system in the nuclear frame {s = = {pi + P2)^), Ni, the incoming 
nucleon and the nucleon in the medium. The factor |piiab|/|pi| is related to the different flux of particles in 
the nuclear frame and in the nucleon frame {pi lab is the incoming LAB momentum in the NN system and pi , the 
momentum in the nuclear system). Furthermore, 



da^ / \P-Po.m.\ \ 



where cm. is the NN center of mass frame and 



0^0 

cm. 



K = xQ{l-\x\) + -^^e{\x\-l), x=^—^^ 



(37) 



(38) 



where Pc.m. is the nucleon momentum in the cm. frame, ei? is the Fermi energy {kp/2M), and the density dependent 
factor Ct is defined in Eq. (A9), with q the momentum transfer in the nuclear frame. In these expressions, 6(k — 
\P-Pc.ni.\) takes into account Pauli blocking and CT(q, p), the nuclear medium renormalization. 



IV. RESULTS 



We compile in Table I the data used for the nuclei studied in this work. As in Ref. [34], nuclear masses and charge 
densities are taken from Refs. [60] and [61], respectively. The neutron density is taken with the same form as the 
charge density but properly normalized and with a different radius as suggested by Hartree-Fock calculations [62] and 
corroborated by pionic atom data [63]. Furthermore, charge densities do not correspond to proton point-like densities 
because of the finite size of the nucleon. This is taken into account by following the procedure outlined in Section II 
of Ref. [63]. 

In the case of NC driven processes, the minimum energy transfer, (7", needed for a proton (neutron) emission 
reaction corresponds to the proton (neutron) separation energy Qp (Qn)- Our Fermi gas model does not account 
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NucIbus 


Rp Rn 


dp dj-^ 




\AZ ) 




D ( A A 




o 


16q 


1.833 1.815 


1.544 1.529 


12.127 


15.663 


- 0.536 


2.489 


14.906 


10.931 


«Ar 


3.47 3.64 


0.569 0.569 


12.528 


9.869 


7.582 


5.830 


0.994 


7.991 




3.51 3.43 


0.563 0.563 


8.328 


15.641 


0.539 


7.799 


13.809 


1.822 




6.624 6.890 0.549 0.549 


8.008 


9.001 


3.707 


3.790 


2.368 


5.512 



(*) The parameter a is dimensionless for the MHO density form. 



TABLE I: Charge {Rp,ap), neutron matter (Rn,a„) density parameters (in fcrmi units), and Qp(„) proton (neutron) energy 
separation, Q— and Q— values for different nuclei in MeV units. For the oxygen we use a modified harmonic oscillator 
(MHO) density, p{r) = po{l + a(r/7?)^) exp(— (r/i?)^), while for the rest of the nuclei, a two-parameter Fermi distribution, 
p(r) = po/(l + exp((r — R)/a)), was used. 

properly for this minimum energy and we correct it by replacing 

9° - 5° - Qp{Az) (39) 

for the proton emission reaction and similarly for the neutron one. In the case of CC processes, this minimum excitation 
energy is Q + Qp{Az+i) [Q + Qn{Az-i)] for neutrino [antineutrino] induced reactions, with Q = M{Az+i) — M{Az) 
[Q = M{Az-i)-M{Az)]. 

In what follows, we will present different results for NC induced inclusive and NC and CC miclcon emission reactions 
at moderate energy transfers to the nucleus, which illustrate the important role played by the different nuclear effects 
considered in this work: Pauli blocking, RPA and FSI effects, and finally the re-scattering of the outgoing nucleon. For 
consistency we will always use non-relativistic kinematics to evaluate the contribution of a particle-hole excitation. 
Relativistic effects were studied in Ref. [34] and found there to be small at the energy regime studied here. 

All differential cross sections shown in this section are computed in the LAB frame. 

A. Inclusive NC scattering at low energies 

In this section we present results for the inclusive QE NC cross sections for both neutrino and antineutrino beams 
in several nuclei concentrating specially on the nuclear medium effects. The results obtained with the same model for 
CC processes can be found in Ref. [34]. 

In Fig. 3 we show the cross section for the processes i'{9) + A — > +X at low and intermediate energies. Above 
the shown range, pion production could become relevant. We also show for comparison the isospin averaged neutrino 
(antineutrino) free nucleon cross section. 

We find that nuclear effects produce a strong reduction of the cross section. The simple consideration of Pauli 
blocking is already quite effective even at the higher studied energies. The inclusion of the long range RPA correlations 
is also quite significant. On the other hand, the magnitude of the reduction is only weakly dependent on the isospin 
and the atomic mass showing, as expected, a larger effect for heavier nuclei. The size of the effect is similar to the 
one found in Ref. [34] for CC processes although in this case we also explore the isoscalar pieces of the effective ph-ph 
interaction responsible of the RPA correlations and the energy thresholds are different. 

The addition of FSI (Sect. II C 2) to the calculation containing already Pauli blocking and RPA correlations does 
not affect practically the results for the totally integrated cross sections. 

We show in Fig. 4 the final neutrino energy spectrum for two typical cases. Very similar results are obtained for 
other nuclei and energies. Although this cross section is not experimentally observable it is useful to show how nuclear 
effects modify the spectrum and whether they favor some energy transfer. We find that RPA produces a quite smooth 
reduction that covers all the energy range, although the effects produced by the RPA correlations are weaker when 
the energy transferred to the nucleus is larger, much the same as it is the case for CC processes. The inclusion of FSI 
spreads the spectrum allowing for larger energy transfers. 

In Fig. 5, we show the ratio of the cross sections induced by neutrino and antineutrinos in oxygen and lead from our 
full model and compare to the results obtained in the free case: LDT (Eq. (19)). The effect of the nuclear medium 
in the ratio is quite small although it increases with energy and with the nuclear mass. Similar results were found in 
Ref. [34] for the CC case. 



Indeed, in a Fermi model we need zero excitation energy to emit nucleons. 
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FIG. 3: (color online). Cross section for the viu) + A p{p) + X processes as a function of the beam energy. Crosses: 
((T„(i7)+p_„(i7)+p + <^v{9)+n^v{u)+n) /"i- Intermediate band: Cross sections for several nuclei including Pauli blocking. Lower 
band: Cross sections for several nuclei with Pauli blocking + RPA correlations. The two solid lines of the lower band correspond 
to the full model calculation (including also FSI) for oxygen and lead. The left and right panel correspond to neutrino and 
antineutrino beam respectively. 




FIG. 4: (color online). NC neutrino da/dE^i cross section as a function of the energy transfer and for two different incoming 
neutrino energies. Predictions from different stages of refinements of the model are shown. Long dashed line: Pauli blocking. 
Short dashed line: Pauli blocking + RPA. Dotted line: Pauli + RPA + FSI considering only the real part of the nucleon 
self-energy as explained in the text. Solid line: full model. 



B. Inclusive {vi,viN), {1^1,1 N), [ui,viN) and {1/1,1^ N) reactions in nuclei at low energies 

In this section we present results for processes, in which the final nucleons are detected for both neutrino and 
antineutrino beams in several nuclei. We wiU also include CC reactions using for this case the model of Ref. [34] in 
which we have have implemented the nucleons rescattering as described in Sec. III. 

As it is illustrated in Fig. 4, the consideration of the dressing of the nucleon propagators (FSI) produces a quenching 
of the QE peak respect to the simple ph excitation calculation and a spreading of the strength, or widening of the 
peak. Most of the effect of the FSI comes from the consideration of the real part of the nucleon selfenergy of Eq. (28) 
as it is shown in the results obtained taking the limit of vanishing ImE (dotted curve in Fig. 4). Thus, we find that 
the change in the nucleon dispersion relation is more important than the inclusion of the small nucleon width in the 
medium, related to the quasielastic channels, which will account for Z'^ absorption by two nucleons. Given the quality 
of this approximation for the production of the spectra needed to start the MC simulation of the nucleons rescattering, 
and the considerable reduction in time for its evaluation, it has been used for all the MC results presented in this 
work. 
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FIG. 5: (color online), a^/ai^ ratio for the viv) + A ^ i^'i^') + ^ reaction as a function of the energy for several nuclei. LDT 
means free nucleon results. In all other lines the full model is used. 



1. CC Nucleon spectra 

The nucleon spectra produced by CC processes induced by muon neutrinos and antineutrinos of 500 MeV are 
shown in Fig. 6 for argon. Of course neutrinos only interact via CC with neutrons and would emit protons, but 
these primary protons interact strongly with the medium and collide with other nucleons which are also ejected. As 
a consequence there is a reduction of the flux of high energy protons but a large number of secondary nucleons, many 
of them neutrons, of lower energies appear. We should recall that our cascade model does not include the collisions 
of nucleons with kinetic energies below 30 MeV. Thus, the results at those low energies are not meaningful and are 
shown for illustrative purposes only in Fig. 6. 

The flux reduction due to the quasielastic NN interaction can be easily accommodated in optical potential calcu- 
lations. However in those calculations the nucleons that interact are lost when in the physical process they simply 
come off the nucleus with a different energy and angle, and may be charge, and they must be taken into account. 

2. NC Nucleon spectra 

The energy distributions of nucleons emitted after a NC interaction are shown in Figs. 7 and 8. In Fig. 7, we show 
the results for '*°Ar at two different energies. In both cases we find the large effect of the rescattering of the nucleons. 
For 500 MeV neutrinos the rescattering of the outgoing nucleon produces a depletion of the higher energies side of 
the spectrum, but the scattered nucleons clearly enhance the low energies region. For lower neutrino energies, most 
of the nucleons coming from nucleon nucleon collisions would show up at energies below the 30 MeV cut. 

As expected, the rescattering effect is smaller in lighter nuclei as can be seen in Fig. 8 for oxygen. In all cases the 
final spectra of protons and neutrons are very similar. Our results without rescattering can be compared with other 
calculations like those of Refs. [22, 32] and other like Refs. [31, 35]. However, in these latter cases, which incorporate 
the nucleon final state collisions, via the use of optical potentials the main effect of rescattering is to reduce the cross 
section at all energies instead of displacing the strength towards lower energies as we find. 

Also interesting is the very similar spectrum shape obtained in shell model [32] and Fermi gas calculations, as the 
present one. 

Finally, in Fig. 9, we show the effect over the spectrum of the RPA correlations and the consideration of the 
real part of the nucleon selfenergy (FSI). We find that the larger reduction due to RPA takes place for the lower 
energy nucleons. A similar situation is obtained for higher energy neutrinos. The inclusion of FSI effects produces an 
enhancement of the number of high energy protons due to the different nucleon dispersion relation. 
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FIG. 6: (color online). Charged current '*''j4r(z/^, ^~ +A'') (upper panels) and '*''ylr(z/^, +A'^) (lower panels) cross sections as 
a function of the kinetic energy of the final nucleon for an incoming neutrino or antineutrino energy of 500 MeV. Left and right 
panels correspond to the emission of protons and neutrons respectively. The dashed histogram shows results without nucleon 
rescattering and the solid one the full model. 



3. Nucleon angular distributions 

The angular distribution of nucleons is also affected by the rescattering. In Fig. 10, we show the proton and 
neutron spectra in the +^*' Ar ix~ + X reaction. For comparison we also show the distribution of protons 
without rescattering. After taking into account the rescattering the proton cross section is less forward peaked. Even 
flatter is the neutron cross section, because neutrons come only from secondary NN collisions and not from the weak 
neutrino-nucleon interaction. 

For NC reactions (Fig. 11), the situation is more symmetric and the angular distribution is similar for protons and 
neutrons, as was the case for the energy spectra. 



4- Energy- Angle distributions 



In Fig. 12-15 we show double differential energy-angle cross sections for both charged and neutral currents. Three 
of the four panels show a common feature with previous calculations (see i.e. Fig. 1 of Ref. [22] and Fig. 5 of 
Ref. [32]), namely, at forward angles there are two peaks, one of them at low energies, that merge into a single one for 
larger values of the angle. As it is shown in Ref. [32] the use of the nucleon momentum distributions from shell model 
wave functions produces less sharp features at forward angles than the Fermi gas calculation of Ref. [22]. However, 
our results are even softer than those of Ref. [32] although our calculation starts with a local Fermi gas momentum 
distribution. The reason is clear, the strong effects of rescattering change direction and energy of the nucleons. This 
can be seen in the panel showing (Fig. 13) the neutron emission for the CC process v^+^^O^ iJ,~ -\-n-\- X . In this case 
all neutrons come from nucleon-nucleon collisions and spread over the available phase space without any remarkable 
feature except the accumulation at low energies. 
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FIG. 7: (color online). Neutral current Ar{u,u + N) at 500 MeV (upper panels) and 150 MeV (lower panels) cross sections as 
a function of the kinetic energy of the final nucleon. Left and right panels correspond to the emission of protons and neutrons 
respectively. The dashed histogram shows results without rescattering and the solid one the full model. 



5. Proton to neutron ratios 

The ratio of proton to neutron QE cross section could be very sensitive to the strange quark axial form factor of 
the nucleon, and thus to the g\ parameter [19, 20, 22, 26, 32]. Our results for this ratio in ^^O are shown in Fig. 16. 
We do not consider very low energies where shell model effects could be more important and our MC simulation is 
unreliable. We find similar results for light nuclei and low energies as in Ref. [20] where RPA correlations were taken 
into account or in Ref. [22] . Our model includes as the main additional ingredient the rescattering of the nucleons via a 
MC simulation. However also this rescattering produces minor changes for light nuclei, because of the smaller average 
density, and for low energies because most secondary nucleons are below our 30 MeV cut. However, the sensitivity to 
the collisions of the final nucleons is larger for both heavier nuclei and for larger energies of the neutrinos as shown 
in Fig. 17 where it is clear that one sees the importance of the secondary nucleons at the low energies side of the 
spectrum. 

V. CONCLUSIONS 

We have studied the QE contribution to the inclusive [ui,uiN), (yi,l~ N), {Di,DiN) and (DiJ^N) reactions in nuclei 
using a MC simulation method to account for the collisions of the ejected nucleons during their way out of the nucleus. 
As input, we have used the reaction probability from the microscopical many body framework developed in Ref. [34] 
for CC induced reactions, while for NC we use results from a natural extension, performed also in this work, of the 
model described in that reference. Limitations of the DWIA models have been discussed. In particular those models 
cannot properly describe individual inclusive neutron and proton spectra and for the ratio of proton (i>, p) to neutron 
(j/, n) yields, the sensitivity to the collisions of the final nucleons might become important for both medium and heavy 
nuclei and for energies of the neutrinos larger than 150 MeV, as shown in Fig. 17. 
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APPENDIX A: NC NUCLEON TENSOR 



1. Impulse Approximation 



Taking into account that in Eq. (17) both the particle and the hole nucleons are on the mass shell {p^ — (p + q)'^ — 
M^, 2p-q + q^ = 0), one finds 

A%^ip, q) = af 5^-^ + + P'^" +P'''i' ^ + ia^e^''" V^/s + afg^q" (Al) 

with, omitting the obvious subindex N = n oi p, 



a,{q') = 8qUiFf+f,zFir + iGir 



4 g2 

asiq^) - WG^iFf + ^IzFi) 

= + i) - ^(^Pff^zFi (A2) 
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FIG. 9: (color online). da/dT^ cross section for oxygen at — 150 MeV as a function of the kinetic energy of the final proton 
for different nuclear models, including or not RPA and FSI effects. In all cases a MC simulation is performed to compute the 
rescattering of the outgoing nucleons. 
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FIG. 10: (color online). Nucleon angular distribution for the Ariv^,, + A*') reaction. The dashed line shows the proton 
results without rescattering and the solid (dotted) histogram stands for proton (neutron) results from our full model. In all 
cases, the contribution of low energy nucleons (below 30 MeV) is not included. 



The cross section for the process vi + N ^ vi + N is given by 

where s — (2|fc | -\-M)M is the Mandelstam variable (|fc | is incoming neutrino energy in the LAB frame). The variable 
is related to the outgoing neutrino LAB polar angle {9') by = {k ~ k')^ — —2\k\\k' |(1 — cos6''). 
The cross section for the process Di + N ^ Di + N is obtained by replacing by —03. 
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FIG. 11: (color online). Proton and neutron angular distribution for the ^^0{v, v' + A'') reaction. In all cases, the contribution 
of low energy nucleons (below 30 MeV) is not included. 




2. RPA Corrections 



Taking q in the z direction and after performing the RPA sum of Fig. 2, we find, neglecting^^ corrections of order 

^ = 8 (^)' {(Cn - 1) mr - (f^f + (Dn - 1) [(Ff r + {F,'n" } 

- 4^{(Cn - 1) [{F^Y - (Ff )"] [(/izi^f - {^izFir] 

+ (Dn - 1) [{F^Y + [FtY] [if^zFiY + i^^zFiY] } (A4) 
^ - (^^^^){4(CN-l)[(^^fr-(Ff)1%4(DN-l)[(Fff + (^^^^^^ 



Note that q^/AI is of the order |gp/M^ and as mentioned in Sect. IIC, we have considered /i^F^lg |/M of order 0(0). 
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RPA 

RPA 
2M2 



(Cl - 1) {{G\r - (Gl)"]' + (En - 1) {{G\Y + (Gl)'f } (A5) 
2(Cl - 1) {[Glf - (Gl)'f + 2(En - 1) {(G\r + (GI)"]' (A6) 

-2^{(Ct - 1) {{^izFir - {i^zFi)T\^ + (En - 1) [(a^zF^^))^ + (a^zF^^))"]' } 

- 4-^{(Ct - 1) {{Fir - (Ff )"] [(mzF^^)^ - {^zFiT\ 
+ (En - 1) {{F^r + (Ff )"] [(mzF/)'' + (/izF/)"] } 

+ 2(Ct - 1) {{Glr - (Gi)"] ' + 2(En - 1) {{Glr + (Gl)"] ' (A7) 

^ _4iM^{(Cx - 1) [(Gl)^ - (Gl)"] [(Ff + MzF/)^ - (Ff + /.^F/)"] 
+ (En - 1) {{G\r + (GDI [(^f + i^zFir + (Ff + MzF#)"] } (A8) 
with the polarization coefficients defined as 

Cn(p) = \rr < , Ms. ^ Ct(p) = „ ] ,,,, Cl(p) = ^ 



|l-co/^(p){/^(g,fcj.)P' |l-[/(g,fcj.)T/*(q)|2 |l-C/(g,fci.)l^z(q)|2 
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FIG. 16: (color online). Ratio of da/dT for protons over that for neutrons for — 150 MeV in the reaction ly+^'^O i/' +N+X 
as a function of the nucleon kinetic energy. Dashed histogram: without nucleon rescattering. Solid histogram: full model. 



Dn(p) = 



1 



\l~Cofoip)UNiq,kFW 



En(p) = 



1 



1 - cogo{p)UN{q, fcp-)!^ 



(A9) 
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